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ble asymmetric nuclear matter even in some cases for "which the symmetric nuclear matter limit 
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hole interaction producing alternative behaviors at different p and b (respectively the density and 
the asymmetry coefficient). The spin and spin-isospin coefficients, with corrections to the usual 
Landau Migdal parameters, indicate the possibility of occurring instabilities with common features 
depending on the nuclear density and n-p asymmetry. Possible relevance for high energy heavy 
ions collisions and astrophysical objects is discussed. 
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1 Introduction 



The symmetry energy coefficients and their dependences on the density and isospin symmetry are 
of relevance, for example, for the description of macroscopic nuclear properties as well as for proto- 
neutron and neutron stars. In the case of the neutron-proton (n-p) symmetry energy coefficient 
(s.e.c) of nuclear matter (a,-) it represents the tendency of nuclear forces to have greater binding 
energies (E/A) for symmetric systems - equal number of protons and neutrons. It contributes as 
a coefficient for the squared neutron-proton asymmetry in usual macroscopic mass formula in the 
parabolic approximation: 

E/A = HQ{A,Z)+ar{N - Zf/A^, (1) 

where Hq does not depend on the asymmetry, Z, N and A are the proton, neutron and mass numbers 
respectively. Higher orders effects of the asymmetry (proportional to (A^ — Z)" for n 7^ 2 are 
usually expected to be less relevant for the equation of state (EOS) of nuclear matter based on such 
parameterizations ||2|, |^. a,- is also the parameter which measures the response of the system to a 
perturbation which tends to separate protons from neutrons. It is given by the static polar izability 
of the system (the inverse isospin screening function) which also should depend on the asymmetry 
of the medium. In the scalar channel one can define the (dipolar) incompressibility which is related 
to the nuclear matter incompressibility [Q]. These two coefficients (isovector and scalar) and their 
dependences on the neutron-proton asymmetry were initially investigated in a previous work [Q]. Other 
symmetry coefficients may also be defined in nuclear matter, for instance, the spin and spin-isovector 
Octt ones. The former would correspond to the difference in the binding energy due the inclusion of 
a polarized nucleon in the medium whereas the second involves the distinction between neutrons 
and protons as well. A calculation for the dynamical and static polarizabilities - proportional to the 
inverse of such symmetry coefficients in asymmetric matter - was done using Skyrme effective forces in 
^ . The spin channel is relevant for the study of the neutrino interaction with matter because it 
couples with axial vector current together with the scalar channel |8|, |^ . A suppression of the spin 
susceptibility (in this work we will be dealing rather with its inverse, the spin symmetry coefficient 
- Aifi, as shown in sections 2 and 3, corresponding also to an increase of the Landau parameter 
Go) lead to the suppression of Gamow Teller transitions which are of interest for the supernovae 
mechanism |^]. Calculations with Skyrme interactions usually result in smaller Go (corresponding to 
smaller spin symmetry coefficients) than microscopic calculations leading to instabilities associated 
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to ferromagnetic polarized states [10, In the present work we show this conclusion about Skyrme 
interactions is not necessarily correct. The spin-isospin channel has been associated with instabilities 
which would lead to pion condensation. In the language of Migdal-Landau parameters one would have 
the parameter Gg < — 1 for the formation of such condensates. 

In the present work we show that these coefficients may provide us with a way of checking the 
stability of asymmetric nuclear matter with respect to the explicit isospin asymmetry. Another aim is 
to articulate the idea that different Skyrme-type forces may be appropriate to the description of diverse 
phenomena at variable ranges of the nuclear density and asymmetry by extending the previous studies 
of s.e.c. Some Skyrme interactions will be used to assess the possible behavior of these functions. The 
parameters of one of the forces we use (SLyb) were fitted from results of neutron matter properties 
obtained from microscopic calculations in Other forces (SkSC4, SkSC6 and SkSClO), which have 
slightly different density dependencies, had their parameters fixed by adjusting a large amount of 
nuclear masses [0, |l^ . 

This paper reviews and extends the works presented in [Q, |^ and it is organized as follows. In the 
next section we remind an argument previously developed to investigate generalized static polariz- 
abities which allow for the study of the nuclear matter stability. In section 3 the static polarizabilities 
in the four channels of the particle-hole interaction with several Skyrme forces is studied in asym- 
metric nuclear matter at variable densities as derived from the linear response method. Next a new 
stability condition with respect to neutron-proton fluctuations asymmetry is proposed. In the last 
section possible consequences are pointed out and the results are summarized. 

2 Generalized Symmetry Energy Coefficients 

We review a qualitative argument from |14, ^ for exploring them. We consider a small amplitude 
(e) external perturbation which acts, through the third Pauli isospin matrix t^, in nuclear matter 
separating nucleons with isospin up and down 0. This originates fluctuations dp = 6pn — Spp = (3 oi 
the nucleon densities. The energy of the system can be written as: 

H = H, + ^0,1 " + 6/3, (2) 



^This argument is valid for the four channels. It is enough to consider other external perturbations: 0-3, 0-3x3 and 1 
for the spin, spin-isospin and scalar channels respectively. 
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where As^t is the isovector symmetry coefficient {s = 0,t = 1 - spin, isospin) usuahy denoted by a,-, J 
or p. For the other channels (s,t) one defines different symmetry coefficients. We can re-write the n-p 
density difference in terms of f3 as: Pn — Pp = iPn,{0) ~ Pp,{0) + P)^ where the subscript (g) indicates 
static densities, without the external source. Although one is dealing with an infinite system, the 
arguments are to be valid for finite nuclei. 

The polarizability is defined as the ratio of the density fluctuation to the amplitude of the external 
perturbation and it can be written as @, U] : 

where we consider the general channel (s,t) as done in j^. Note that As,t is a function of b and /? and 
these parameters may be related, as argued below. The occurrence of these functional dependences of 
As^t can be found just by the stability condition with respect to the fluctuation from expression (|2|): 

2.1 Isospin dependence of As,t 

The neutron proton asymmetry used in the present work is defined by the neutron and proton densities 



Pn-i Pp as: 



6 = ^-1. (5) 

Pp 



It varies from 6 = 0, in symmetric nuclear matter, up to infinity, in neutron matter. The coefficient b 
is related to a frequently used asymmetry coefficient: 

a = 1^^211^, (6) 
Po 

by the expression: 

b = 2a/(l - a). (7) 

As stated above As,t is a function of the density fluctuation /3. Although (3 is not the explicit n-p 
asymmetry itself (given by b) we will consider that it depends on it (as it was also argued in [Q, ^]). 
We consider, as shown below, these parameters are related to each other and therefore we will write 
A = A{I3) shortly. In |5|, |6| different prescriptions were discussed for (3 = (3{b) in the calculation of the 
response function of asymmetric nuclear matter. We have used (and it was shown to be a reasonable 
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prescription for the dynamical response function the one which leads to the following relation 
between the fluctuation (3 and the explicit n-p asymmetry h: 

/3 = <5Pnf^V (8) 



,1 + 6 

Where 5pn is the neutron density fluctuation. In the n-p symmetric limit (3 = 2dpn and in another 
limit, in neutron matter, /3 = 6pn- The above prescription (expression (|8|)) is based on the assumption 
that the density fluctuations are proportional to the respective density of neutrons and protons, i.e., 
Spn/P = Pn/ Pi being p the total density. In spite of being rather appropriated for the isovector channel, 
this kind of assumption can be considered as a starting point for the other channels (spin, scalar) in 
asymmetric nuclear matter. Prescription (^) is therefore model-dependent and different choices for it 
yield other forms for the (asymmetric) static screening functions. The dynamic response functions are 
less sensitive to this prescription 

^From the solution of the polarizability ^ we calculate the first derivative with relation to h: 

dp _ ep dA^'^ P_cL^ 

lib ~ 2Ali db ~ ~Ao,i db ■ 

Another expression can be obtained from the relation between b and /3 of (|8|). It yields: 

dp _ P 

d6 ~~(2 + 6)(l + 6)' ^^^^ 
Equating these two last equations we obtain: 

(2 + 6)(l + 6) db ' ^ ^ 

^From which it is possible to derive the following relation between the isospin s.e.c. and the n-p 



asymmetry Q: 



In this expression Agym = ar — 30MeV is the s.e.c. of symmetric nuclear matter (6 = 0). For 6 = 2 
(neutron density three times larger than the proton density) we obtain A = 1.5Asym- In the limit of 
neutron matter A{b oo) = 2Asym- The prescription of expression (^) was the relevant information 
for this calculation. Any other relation between 6 and P will induce different asymmetry dependence 



of the symmetry energy coefficient Another assumption for deriving expression (12) was that p 
is independent of 6. This would be non trivially different if one considers a complete self consistent 
calculation with the equation of state of a proto-neutron star, for example. 
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3 Generalized Screening functions with Skyrme forces 

In this section in in the next one we nearly extend the analysis of [^. A nearly exact expression for 
the dynamical polarizability of a non relativistic hot asymmetric nuclear matter at variable densities 
was derived with Skyrme interactions in ||5|. The approximations were: (i) to equate the asymmetry 
coefficient defined for the momentum density to the density asymmetry coefficient (the variation of 
the corresponding term yields differences in the response function of less than 1%), (ii) to choose a 
particular (discussed and justified H]) prescription for the asymmetry density fluctuations - expression 

The generalized static screening functions Ag^t (or symmetry energy coefficients according to the 
discussed before) in asymmetric nuclear matter at finite temperature can be written as ||5|: 



(13) 



{vl^'f [^Q{M;fpcPd - iqm;McN^]. 

The densities pv, My and Ny are given respectively by: 

Pv = vpn + (1 - v)pp, My = vMn + (1 " v)Mp Ny = vN^ + (1 - v)Np, 

where v stands for n-p asymmetry coefficients (c, d) defined below (a measure of the fraction of neutron 
densities). The above densities at finite temperature are defined by: 

2M* r 

{Ng,Pg,Mg) = —f j df ,{k) {k , k\ k^) . 

In these expressions fq{k) are the fermion occupation numbers for neutrons {q = n) and protons {q = p) 
which will be considered only for the zero temperature limit. In this case: dfq{k) = 5{k — kp)dkq, 
which makes the above integrals trivial. Vq and Vi are functions of the Skyrme forces parameters (see 
in |5|) and M* = m*/{l + a/2) is a modified effective mass for the proton. Besides that, the four 
asymmetry coefficients are: 

a = — -1, b= 1, c=——, d = (14) 

ml Pop 2 + 6 1 + (1 + 5)3 

3.1 Results for Skyrme interactions 

In this section we present results for the s.e.c. using Skyrme forces with the expression (|T^ ) in all the 
four channels of the particle-hole interaction. The scalar, isovector and spin channels were initially 
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discussed in a previous works y] and those analysis is extended here. The present analysis on the 
spin and spin-isospin channels was not done before. 

In Figure 1 we show the generalized isovector symmetry energy coefficient Aq^i as a function of the 



ratio of the nuclear density to the saturation density (p/po) for Skyrme interactions SLyb |11], SkSC4, 
SkSC6 and SkSClO [l^. Since Skyrme forces are not necessarily expected to describe Physics at 
high densities we investigate the behavior of the symmetry energy coefficients (s.e.c, which can also 
be called generalized screening functions for each channel) from 0.5 po up to 2 pQ. The symmetry 
energy coefficient for some of the forces in Figure 1 reach a maximum value when po < p < 1.7/>o and 
then decrease. This last behavior is typical of non relativistic calculations. There are two forces which 
yield different behaviors: SkSC4 makes the slope of Aq^i much less negative reaching an instability 
point (^0,1 < 0) at smaller densities than the other forces and SkSClO for which the slope is larger and 
positive. This last behavior is present in relativistic calculations. The point in which the symmetry 
energy coefficient is zero would correspond to a phase transition of nuclear matter to an asymmetric 
state of neutron or proton matter. Although the usual Skyrme forces are not well suited for high 
densities there is a tendency of a decrease of this screening function (s.e.c.) until a phase transition in 
a dense nuclear system. The dependence on b, on the other hand, leads to a more repulsive behavior 
as seen in Figures 2 and 3 eventually compensating the attractiveness of low or high densities. 

In Figure 2 the same parameter, Aq^i, is shown as a function of p/po for the interactions SLy, SkSC4 
and SkSC6 for non zero n-p asymmetries, i.e., for b=.25 and b=.54. The latter corresponds to a n-p 
asymmetry of the nucleus of ^^^Pb. The absolute values of the isovector symmetry energy coefficient 
(s.e.c) increase with relation to the symmetric nuclear matter but the behavior with varying p is nearly 
the same in general. However, for the force SkSC6, the isovector screening function has a higher slope 
in this range of densities. The same coefficient as a function of the n-p asymmetry b is shown in Figure 
3 at the densities 0.5, 1 and 2pQ. Again we note different behaviors depending on the force. For the 
forces SLy and SkSC4 the s.e.c. increases, as noted above, for higher densities and n-p asymmetries 
but for p = 2pQ the values are smaller at higher asymmetries. For the force SkSC6 the higher the 
density the higher the slope with b. The tendency is that the interaction becomes more repulsive. 
Whereas the dependence of Aq^i on 6 is a new possibility its dependence on p has been studied quite 
extensively and a distinction between the behavior of relativistic and non relativistic models can be 
done which still are nearly observed here. While relativistic calculations lead to a continuous increase 



of the symmetry energy coefficient |15|, non relativistic (microscopic or not) calculations usually yield 
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a saturation at densities of the order of 2.5 po with a decreasing value at higher densities |15, 16 1. 
Although these are the more common results there is a non relativistic variational calculation with 
three body forces which presents an increasing symmetry energy coefficient for higher densities. 



All these calculations adopt the parabolic approximation for the symmetry energy as discussed by [17 
i. We have just shown, however, that it is possible to obtain the behavior typical from relativistic 
calculations with the Skyrme force SkSClO. 

The spin symmetry energy coefficient is plotted in Figures 4, 5 and 6 as a function of p/po and b 
respectively. In Figure 4 one sees that, for increasing densities, Ai^ may have different slopes. For the 
interactions SkSC4, SkSC6 and SkSClO the spin- symmetry energy coefficient is already negative at 
low densities and decrease still more with p. This feature is clear in Figure 5 where the same spin s.e.c. 
is shown for non zero n-p asymmetries. These forces exhibit a behavior typical of a ferromagnetic 
matter with stronger spin alignment as density increases. With the force SLy, on the other hand, the 
spin symmetry coefficient indicates that no polarized state occurs the range of densities under study. 

In Figure 6 the spin symmetry energy is shown as a function of b with forces SLy, SkSC4 and SkSC6 
for some different densities. The common trend is the increase of Ai^ with b, i.e., at very asymmetric 
matter the spin interaction tends to become more repulsive. However the particular behavior of the 
spin s.e.c. with b is different for each effective force at a given density. This means that whereas for 
the SLy force the increase of n-p asymmetry makes ^i^o to increase (proportionally) at each density 
considered, for the forces SkSC4 and SkSC6 the slope is higher for higher densities. This makes the 
curves of ^i^o for different densities to cross at certain n-p asymmetry. We can compare our results to 
the ratio of spin susceptibility of interacting neutron matter (for which 6 — > oo) to the non interacting 
Fermi gas obtained by Fantoni, Sarsa and Schmidt ||l^ by means of the auxiliary field diffusion Monte 
Carlo method. This ratio is proportional to the polarizability as obtained in expression ( p!3| ) and 
therefore inversely proportional to the spin symmetry coefficient ^i,o- First of all we note that, in 
most cases, contrarily to what we find, the values they find are all positive for the range of densities 
considered by them, from 0.75po up to 2.5po. However the slope seems to be nearly the same as that 
we obtain for low values of the n-p asymmetry. Consequently they may obtain instabilities for higher 
density neutron matter whereas we do not observe this result in our calculations with Skyrme forces 
(the comparison is meaningful for neutron matter: b is very large). Other comments about this are 



^We leave for a forthcoming work a discussion about corrections to the parabohc approximation in the frame of a 
relativistic model, following ideas contained in section 2 of the present work and in Q . 
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drawn in section 5. 

The spin-isospin symmetry energy coefficient for symmetric matter (b = 0) as a function of p/po 
is shown in Figure 7. For all the forces, Ai^i are decreasing functions of the density and they become 
negative between p ~ 1.8po and 2.1pQ. The only significant difference comes from the value of ^i^i 
at low densities, being that SLy force implies a larger value than the SkSC forces. In Figure 8 the 
same analysis holds for asymmetric matter for b = .25 and .54. However, with increasing asymmetry, 
the instability point is changed differently for each of the forces. This is very interesting as well and 
very similar to the spin channel (of figure 5). The SLy force presents the spin-isospin coefficient which 
changes much less than the SkSC forces. For these Skyrme interactions instabilities associated to pion 
condensation would be found already at low densities, Ai^i < . Therefore, for the Skyrme type 
interactions we use, we expect a there is a connection between the appearance of spin instabilities 
(leading to ferromagnetic states) and spin-isospin instabilities (for pion condensation), even for not 
very asymmetric nuclear matter (up to 6 ~ 0.54). 

In Figure 9, Ai^i is plotted as a function of b for different Skyrme forces. In this Figure the 
behavior of the spin-isospin s.e.c. with the n-p asymmetry parameter b is shown for different nuclear 
matter densities (below, equal to and above the saturation density). It can be noted that for the force 
SLy an increasing n-p asymmetry (b) makes the spin-isovector interaction continuously more repulsive 
instead of allow for attractiveness while the opposite behavior is found for SkSC6. One would find the 
instabilities associated to pion condensation at roughly any nuclear density for SkSC6 at nearly any 
n-p asymmetry, which does not seem to be realistic The other forces SkSC also reproduce this 
last behavior with different slopes. We are lead to conclude that Skyrme interactions may be able to 
reproduce roughly well any kind of dependence of the (generalized) symmetry energy coefficients in a 
particular range of density and/or n-p asymmetry. 

In Figures 10 and 11 the scalar s.e.c. is shown as a function of p/po for b=0 and b=0.25 and 
.54 respectively. These Figures show that, roughly speaking, all the interactions are attractive or 
weakly repulsive at low densities and more repulsive as density increase. This may be related to the 
instabilities found, for example, in [^]. However it is interesting to see in Figure 11 that, depending 
on the interaction, with the increase of the n-p asymmetry the scalar symmetry energy coefficient is 
reduced. This occurs for lower densities but not necessarily at higher p. This is well seen in Figure 
12 for Aq^q as a function of b. This may be understood as a general tendency of resulting instabilities 
for (very) asymmetric nuclear matter with a curious exception for force SLy at a density 2po- The 
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behavior of the incompressibihty with increasing n-p asymmetries, at the saturation density, for the 
force SLy is nearly the same as that found in other calculations with Skyrme forces or relativistic 
calculation ^ and microscopic variational study at finite temperatures |^^. For a higher density 
(2/9o)) however, this scalar coefficient increases with b for the same SLy force. A comparison of this 
kind was explicitly done in Q. The other forces exhibit smaller slopes and make Aq^ to decrease 
faster. This last kind of behavior, more pronounced, was also found in other works (for the bulk 



incompressibihty modulus K^o depending on the used interaction [22, Therefore the behavior of 
the incompressibilities (i^oo and Aq^ - the "dipolar incompressibihty") with n-p asymmetry depend 
strongly on the used interaction. 

4 Stability of Asymmetric nuclear matter 

Consider that the binding energy can be minimized with respect to /3 yielding an equilibrium state 
for nuclear matter. This means that: d(E/A)/dp = 0, yielding a differential equation for As,t as a 
function of /5 or Hs^t by expression (^). In the case p is not dependent on /5 the solution for this 
resulting differential equation is given by Q: 

As,t _ C 1 

~^ ~ TT2 ff ' 

where C is a constant. Inverting the above expression we find: 



2A,^^ As, +4^2^' 

which, in principle, can be either real or complex. Furthermore, there may appear two branches 
of solutions for 6 7^ 0. We want to emphasize that the relation between this expression and those 
obtaining from the linear response with Skyrme forces is not completely understood. As a boundary 
condition one requires the usual expression for the symmetric limit obtaining: 

AA7^r' ^ ' 

Such that the corresponding polarizability of symmetric nuclear matter IIq* by: IIq* = —pq/{2A7JI^) 



To really be a stable minimum of the binding energy the second derivative of the energy must be 
positive. Using the form given by expression (^), and the approximation that pQ does not depend on 
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P, we can write this second derivative as: 



d^H 1 4UdA U^d'^A ,,,, 
^ + ^ + (18) 



d/?2 U po dU po dn2 
With the solution given by expression ([l5| ) we find, in the (s, t) channel, that: 



The constant C may be negative (stable symmetric nuclear matter according to expression (p!7|)) or 
positive and A and 11 also may have negative or positive signs. It is possible to have a case where the 
symmetric nuclear matter is unstable, Agym < 0, but a stable asymmetric nuclear matter is obtained, 
with positive d^H/dP^. This is a curious result. It seems to suggest that the stability line of nuclei 
can occur, at least with its qualitative real behavior, without introducing explicitly electromagnetic 
forces. 

For the sake of comparison we quote the usual stability condition for a Fermi liquid in each channel 
of the interaction. It is given by: 

as,t = iVo(l + Jo'*) > 0. (20) 

where Jg'* stands for any of Fq, Fq, Gq, G'q respectively for the scalar (s = 0, t = 0), isovector (s = 0, t = 
1), spin (s = l,t = 0) and spin-isovector (s = l,t = 1) channels |2^. These conditions correspond to 
the denominators of the response function of symmetric nuclear matter at zero temperature and at 
saturation density. 



5 Further discussion and Summary 

The expressions we have derived and used for the (isospin) symmetry energy coefficient (with and 
without effective Skyrme forces) break isospin symmetry because the results for an excess of neutron 
(density) is different from those of higher proton fraction. In principle this is expected although 
usually considered to be small in mass formulae for finite nuclei [Q]. With the expressions we derived 
the effects of this symmetry breaking can be expected to be larger eventually related to other effects 

Recently it has been proposed that in the core of dense (neutron) stars there may occur a phase of 
matter in which color degrees of freedom are deconfined generating di-quark condensates in channels 
in which the interaction is attractive. Several possibel scenarios for this color superconducting phase 
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of matter can be drawn, among them the color-flavor locked phase [^]. We would like to point out 
that with the increase of the isospin symmetry energy coefficient, for instance, the resulting fraction 
of protons, after the supernovae stage, would be higher hindering the emergency of a pure neutron 
matter in the death of the star. As a consequence the final proportion of up and down quarks would 
be different than that from a pure neutron star. This may modify the color superconducting phase. 
These remarks may apply to relativistic or high energy heavy ions collisions which generates high 
baryonic density. 

The relevance of the isovector symmetry energy coefficient and its dependence on the nuclear 
density has been intensively studied due to its relevance to the determination of the nuclear matter 
equation of state which is very important for the comprehension of high energy ions collisions and 
eventually the dynamics of astrophysical objects. An interesting recent study was done by Bao An 
Lin who however only is concerned about the parabolic approximation without a more general 
form of the asymmetry term. This extension would be of high interest for the field. 

Summarizing, the dependence of the s.e.c. on the n-p asymmetry was studied extending the 
results of ref. |Q, ^ . Generalized symmetry energy coefficients of asymmetric (non relativistic) nuclear 
matter at variable densities were investigated. The density and n-p asymmetry dependences of the 
s.e.c. in the different channels were analyzed for different Skyrme forces. They may yield very 
different behaviors including the possibility (or not) of nuclear matter to undergo phase transitions. 
These forces can describe different behaviors of the symmetry energy coefficients. Although one 
should not believe that only one Skyrme force parametrization could account for the description of all 
nuclear observables at different densities, temperatures and n-p asymmetries it is acceptable the idea 
that several parametrizations could hopefully describe different ranges of the dependence of nuclear 
observables (as the s.e.c.) with these variables. Whereas the isovector, scalar and spin channels were 
rather an extension of the work presented in [Q, |6| the results for spin and spin-isospin channels are 
scarce in the literature. We showed that, for the Skyrme forces under consideration (SLy, SkSC4, 6, 10) 
there is a kind of correlation between the two channels and the possibility of occurring simultaneous 
ferromagnetic instabilities and/or pion condensation at high densities and n-p asymmetries. Therefore, 
in principle, different values can be expected for the (bulk) symmetry energy coefficients in nuclei and 
nuclear matter, with different n-p asymmetries, and neutron stars. In the spin channels it is possible 
to expect spin polarized asymmetric matter yielding magnetic fields in neutron stars, as discussed in 



|10| , according to the results in figures 4,5 and 6 for forces the SkSC. However with the increase of 
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b we find that the spin interaction may be rather repulsive, hindering this magnetization effect with 
the use of these Skyrme forces in the range of densities analyzed here. A similar behavior for the 
spin-isospin channel (which may indicate instabilities for pion condensation) was found with the same 
Skyrme interactions. A new stability condition for an asymmetric nuclear medium was proposed. It 
is given by: 



Where C is a constant and 11*'* = p/ {2As,t) for each channel. 
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Figure 1 Neutron-proton symmetry energy coefficient ^0,1 = p/(2n^ ) of symmetric nuclear matter 
as a function of the ratio of density to the saturation density (n = p/po) ^01 interactions SLyb (dotted- 
dashed hne), SkSC4 (solid), SkSC6 (dotted), SkSClO (dashed). 

Figure 2 The same as Figure 1 for interactions SLy with b=0.25 (thick dotted-dashed line), b=.54 
(thin dotted-dashed line), SkSC6 with b=.25 (thick dotted), b=.54 (thin dotted), SkSC4 with b=.25 
(thick solid) and b=.54 (thin solid). 

Figure 3 Neutron-proton symmetry energy coefficient ^0,1 = p/ (211^^) as a function of the asymmetry 
coefficient b at different densities: solid lines for SkSC4 (very thick: p = .5po, thick:/9 = pQ, thin:p = 
2pQ,), dotted lines for SkSC6 (same conventions as SkSC4) and dotted-dashed lines for SLy. 
Figure 4 Spin symmetry energy coefficient Ai^ = p/{2Il^^^) as a function of the ratio of density 
to density at saturation {u = p/po) for interactions SLy, SkSC4, SkSC6 and SkSClO with the same 
conventions of Figure 1. 

Figure 5 Spin symmetry energy coefficient Ai^ = p/{2Il]f) as a function of the ratio of density 
to density at saturation {u = p/po) for interactions SLy, SkSC4, SkSC6 and SkSClO with the same 
conventions of Figure 2. 

Figure 6 Spin symmetry energy coefficient Ai^q = p/ {211^^) as a function of the asymmetry coefficient 
b at different densities and for the different forces: with the same conventions of Figure 3. 
Figure 7 Spin-isospin symmetry energy coefficient Ai^i = p/(2n^^) as a function of the ratio of 
density to density at saturation {u = p/po) for interactions SLy, SkSC4, SkSC6 and SkSClO with the 
same conventions of Figure 1. 

Figure 8 Spin-isospin symmetry energy coefficient Ai^i = p/(2n^^) as a function of the ratio of 
density to density at saturation (n = p/po) for the different interactions: with the same conventions 
of Figure 2. 

Figure 9 Spin-isospin symmetry energy coefficient Ai^i = p/{2Il]^^) as a function of the asymmetry 
coefficient b at different densities and with the different forces: using the same conventions of Figure 
3. 

Figure 10 Scalar symmetry energy coefficient Aq^ = p/{2Il^) as a function of the ratio of density 
to density at saturation (n = p/po) for the different interactions with the conventions of Figure 1. 
Figure 11 Scalar symmetry energy coefficient ^0,0 = p/(2n^°) as a function of the ratio of density 
to density at saturation {u = p/po) for the different interactions: with the conventions of Figure 2. 
Figure 12 Scalar symmetry energy coefficient Aq^ = p/{2Il^jf) as a function of the asymmetry 
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coefficient b at different densities and with the different forces, the conventions of Figure 3. 
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Figure 4 
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Figure 7 
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Figure 10 
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